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It is shown using three series of Rayleigh number sim- 
ulations of varying aspect ratio AR and Prandtl number Pr 
that the normalized dissipation at the wall, while significantly 
greater than 1, approaches a constant dependent upon AR 
and Pr. It is also found that the peak velocity, not the mean 
square velocity, obeys the experimental scaling of Ra 05 . The 
scaling of the mean square velocity is closer to Ra 46 , which 
is shown to be consistent with experimental measurements 
and the numerical results for the scaling of Nu and the tem- 
perature if there are strong correlations between the velocity 
and temperature. 



This Letter will analyze the energy budget in three- 
dimensional simulations of Rayleigh-Benard convection 
with the objective of testing theoretical assumptions used 
to explain the laboratory observation of non-classical 
heat flux (Nusselt number Nu) exponents in classical tur- 
bulent Rayleigh-Benard KM- The budget equation for 
the average kinetic energy <p /2 as a function of height z 
is 



1 d^r 

2di q = 



-pw z — 0.5wq 2 z + RawO + Prd z q 2 — Pre(z) 

(1) 



where shear production 



uw(dU/dz) is part of the tur- 
Based upon_the ob- 
and 



bulent production term — O.bwq 2 z 
servation of strong shears in the boundary layer 
large-scale flows B, it was first assumed that the dissi- 
pation is primarily concentrated in the boundary layer 
and is turbulent M . More recently it has been suggested 
that the boundary layer is laminar and the distribution 
of the total energy dissipation is Rayleigh number depen- 
dent [§J. Both of these approaches predict crossovers in 
scaling behavior. Another theory makes mixing layer 
assumptions and does not predict these crossovers. The 
objective in this Letter is to look at the basis for these 
assumptions using numerical simulations. 

Assuming that Nu — wO / ' ndQ / ' dz ~ RaP T , the clas- 
sical exponent is (3r = 1/3, with suggestions since the 
mid-60 's 12 R| that there might be corrections to this 
exponent. The original experimental result [jj] showing 
that there are significant corrections has recently been 
extended and refined to give Nu ~ Ra 309 over nearly 10 
decades of Rayleigh number Ra However, detailed 
experimental information such as budgets cannot be ob- 
tained from these large Ra experiments. The only in- 
formation available besides temperature statistics at a 
single point is that the Reynolds number based upon 



vertical velocity fluctuations w midway up a sidewall 
goes as Re = wd/v ~ Ra 1 / 2 . These two observations 
are incompatible with standard turbulent parameteriza- 
tions because the total dissipation is constrained to be 
€t = Ra(Nu — 1) ~ Ra 131 , while the standard turbulent 
prediction of = Re 3 gives er ~ Ra 15 . 

One reason a laminar boundary layer has been sug- 
gested is that this incompatibility can be resolved if the 
laminar relationship for dissipation is used, ex ~ i?e 5 / 2 ~ 
Ra 1/2S . This would suggest that convective flows are not 
filled with cascading eddies, but are instead filled with 
strong local laminar shears, including a laminar bound- 
ary layer. One way to possibly determine whether the 
boundary layer is laminar or turbulent is to divide the 
total energy dissipation er at some arbitrary boundary 
layer thickness Xbl into the dissipation in the boundary 
layers €bl and the dissipation in the bulk es- A new 
theory for convective scaling M has assumed that in a 
turbulent boundary layer that €bl will be the same or- 
der or greater than es, while if the boundary layer is 
laminar then cbl/cb will decrease as Ra increases. This 
is justified using the scaling laws above for dissipation in 
laminar and turbulent boundary layers. 

Since the average (or total) dimensionless dissipation 
et = Ra(Nu— 1) is increasing with Ra, it is possible that 
even if the boundary layer is laminar at one Ra, that an 
Ra could be reached where the boundary layer becomes 
unstable and turbulent. Then for higher Ra, as in a shear 
driven turbulent boundary layer, cbl/^b — >constant 
would appear. 

One way to obtain detailed information and test these 
properties is to use simulations. Simulations have repro- 
duced most of the scaling laws and statistical prope rties 
of the higher Rayleigh number experiments [|14| , |l5| , in- 
cluding /3t between 2/7 and 1/3. One of the properties 
that simulations can determine are thermal and velocity 
boundary layer thicknesses. The thermal boundary layer 
thickness At can be determined using either l/Nu or 
the position of the peak of temperature fluctuations 9 2 . 
Based upon experience from analysis of classical shear 
driven boundary layers, three definitions of the velocity 
boundary layer thickness Xbl that can be calculated for 
convection use the velocity to give z* = 1/Re, use the 
wall shear stress r = du/dx to give z + = y/r, and use 
the position of the peak of the horizontal velocity fluc- 
tuations X u . The dimensionless forms have been used in 
these definitions. 

These and additional definitions of the boundary layer 
thickness could be used to define the dissipation in the 
boundary layer 6bl, which we would like to relate to 
the dissipation at the wall. While ebl depends upon 
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the definition of the thickness of the boundary layer, a 
relationship that could hold for all the definitions above 
is 

£bl/£t < (ew^BL/(de T )) (2) 

where ew is the dissipation at the wall. This relation- 
ship is not very restrictive, but if cw/ct were bounded, 
because Xbl is a decreasing function of Ra for all the 
definitions given above, €bl/^t would be bounded. 

In a classical shear-driven boundary layer, all of the 
definitions of Xbl given above scale with Re in the same 
manner. A major difference in convection simulations 
fj^flEf is that each definition of Xbl has a separate power 
law dependence on Ra. That is, if 

A u ~ Ra'^ and z+ = ^ ~ RaT^ (3) 

then {3 U w 1/7 and /3+ > /3 T has been found Q. Be- 
cause these results and the theories all indicate that Xbl 
is decreasing rapidly with Ra, by applying (Q) both the 
turbulent and the laminar boundary layer theories pre- 
dict that ew I ct should increase as Ra to a power law. 

Three simulations have been analysed to determine the 
dependence of ew/^T upon Ra. The cases to be dis- 
cussed are AR = 4 and Pr = 0.3 for 10 5 < Ra < 10 7 , 
AR = 4 and Pr = 7 for 10 4 < Ra < 10 7 from earlier 
work p[|, and AR = 1 and Pr = 0.7 for 10 6 < Ra < 
8 x lO^where Prandtl number Pr = v/k and aspect 
ratio AR — width /height. The numerics [ [i"5| are pseu- 
dospectral, using Chebyshev polynomials in the vertical 
to provide more resolution and no-slip, constant temper- 
ature boundary conditions at top and bottom walls. In 
the horizontal, sines and cosines are used to represent 
free-slip, insulating walls. Profiles will be shown only for 
the AR = 1 case to save space and because this case is 
new. 



0.90 



0.85 







-1 


Ra = 8e7 


Pr=.7 






i-i. 




^ i 1 1 


















e N - 














<0w> N 










1 .0 
0.9 




— <wqq/2> 


,z 


; 1 
1 
- 1 












- 1 


i 2 | ,! 


5 10 1- 


5 



-10 1 2 3 

e N <6w> N 



FIG. 1. Production, dissipation and turbulent transport of 
kinetic energy from (hh for the simulation AR = 1, Pr = 0.7 
and Ra = 8 x 10 7 . All dimensionless terms are normalized by 
RaNu so that the production across the center w8n , which is 
also the heat flux, is 1. Insets show production and dissipa- 
tion through the entire box (height d = 2) and very near the 
wall, where dissipation is very large. The difference between 
production and dissipation w6(z) — e(z) is the total transport 
term. As discussed, the pressure transport is close to the 
difference between the production and dissipation and so is 
much larger than the turbulent transport wq 2 /2 . 

Previous work |l4f| demonstrated that ew was much 
larger than the average dissipation and that ew / £t was 
increasing with Ra. Fig. [l] shows some of the terms in eq. 
(|Tj) for AR = 1, Pr = 0.7 and Ra = 8 x 10 7 , all normal- 
ized by RaNu. For example, e^{z) = e{z) / RaNu. The 
AR = 4 cases are similar except that ewN — £w I 'RaNu 
is smaller, min(ejv(z)) is closer to 1, and wq 2 /2 is even 
smaller than here. 
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FIG. 2. ewo(AR, Pr) — ejv vs. Ra. The decrease is approx- 
imately Ra~°' s±lhl for all three cases (neglecting the lowest 
Ra for Pr = 0.3 and 7). The values of ewo{AR, Pr) used are 
given. The inset shows the dependence on aspect ratio of Nu 
for Ra = 10 7 and Pr = 0.7. 

While this analysis has confirmed that ew is much 
larger than the average dissipation and ew / £t is increas- 
ing with Ra, careful examination revealed that there 
appears to be an upper bound to ew/^T, which will 
be denoted ewo and is found to be strongly dependent 
upon Pr and AR. Fig. || shows the normalized dissipa- 
tion at the wall plotted as ewo — £wn- For all three 
cases, ewo — € wn ~ Ra~°- 8±01 . While the choice 
of this particular form is subjective, the consistency in 
£\vn — ► ewo (AR, Pr) as Ra grows for all AR and Pr 
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appears to be robust and contradicts the predictions of 
both the turbulent and laminar boundary layer theo- 
ries. Due to the bound (||) on cbl, then cbl/cb — ► 
as Ra — > oo at least as fast as 



£bl/cb < cwoXbl ~ Ra 



-1/7 



(4) 



found p6[ , and confirmed by these simulations, that there 
is a nearly perfect correlation in the bulk between the ver- 
tical velocity w and temperature fluctuations 8, where 
the total temperature is T(x,y,z) — T(z) + 9(x, y, z). 
What is found everywhere in these calculations, except 
near the wall (z < Xt = d/Nu), is that 



where —1/7 comes from using Xbl = X u from (y). 

If this trend were to continue to higher Ra it would im- 
ply that the boundary layer should not be characterized 
as either a laminar or a turbulent shear-driven bound- 
ary layer as has been assumed up to now. Furthermore, 
it implies that the smallest length scales in the prob- 
lem are multiples of the Kolmogorov scale rj — e T ^ 4 = 
(Ra(Nu — 1))~ 1//4 . For example, the wall boundary layer 
thickness taken from the wall shear stress z + = %J~t 
should scale as r). If fa = 0.309, then f3+ = 0.327 is 
predicted. 

Now let us consider the mechanisms responsible for 
transferring energy from the bulk to the boundary layer. 
This is necessary because Fig. [I] shows that the produc- 
tion of kinetic energy, which is equivalent to the convec- 
tive heat flux w0, is found only in the center of the box, 
but the peak of the dissipation is at the wall. For there 
to be a turbulent boundary layer, there would have to 
be large turbulent production or shearproduction terms. 
With the normalization used in Fig. [l| w8/ (RaNu) w 1 
in the center and e(z)/(RaNu) is slightly less than 1, 
which is compensated for by extra dissipation in the 
boundary layer. 

The three transport terms in ([!]) that transfer energy 
from the bulk to the boundary layer are diffusive trans- 
port Prd 2 q 2 , the pressure transport — pw z and the tur- 
bulent transport — 0.5wq 2 z . Prd 2 q 2 and — 0.5wq 2 z can 
be calculated directly while — pw z can be calculated from 
the difference between all the remaining terms. Diffusive 
transport (not shown) is found to be large only very near 
the wall, for l-z < 0.02 in Fig. 0. For l-z < 0.1 the tur- 
bulent transport is much less than the difference between 
the production and the dissipation, and therefore the tur- 
bulent transport is much less than the pressure transport 
term. This would be consistent with the dissipation in 
the boundary layer being much less than either the pre- 
dictions of a turbulent or a laminar boundary layer, both 
of which depend upon shear production, which is part of 
the turbulent production term. 

The observation that the pressure transport dominates 
is not new. In the geophysical literature over a wide 
range in Ra beginning at low Ra in the laboratory ||] 
and extended to higher Ra in atmospheric observations 
0,0, it is found that the primary source of kinetic 
energy in the boundary layer is the pressure transport 
term. 

If the boundary layer is neither a shear driven turbu- 
lent nor a laminar boundary layer, then what mechanism 
is responsible for the observed scaling of the velocity scale 
or Reynolds number Re — wdjvl Perhaps the origin is 
in the details of plume dynamics. It has previously been 
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(5) 



This would be consistent with visualizations that show 
small plumes dominating in the boundary layer J^,[w],[^| . 
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FIG. 3. Reynolds numbers based upon the mean horizontal 
velocity in the direction of maximum shear Vd/v, upon the 
horizontal fluctuations in velocity ud/u — {u 2 + v 2 ) 1 ^ 2 d/u, 
upon the vertical velocity fluctuations wd/v, and the temper- 

The inset shows the profiles 



ature fluctuation variance f? 2 ^ 2 
over the entire domain. 



Using (g), it can be predicted that in the center that 
if Nu = e^/ K (dQ/dz) - Ra t3T , 6' = ¥ ±/2 - Ra' 6 -, and 

—1/2 



Re 



jvd ~ Ra 1 that 

7 = f3 T + S c . 



(0) 



The scaling for the simulations and experiments dis- 
cussed here is given in Table |[ For the experimental 
values § of fa = 0.309 and S c = 0.145, then (g) would 
predict 7 = 0.454, which is the origin of the value given 
in Table 1. However, what is found (private communica- 
tion) for the most recent experiment || and most earlier 
experiments jlTj is 7 = 0.5. For AR — 4, Pr — 0.3 and 
AR = 1, Pr = 0.7: 7 = 0.46 and for AR = 4, Pr = 7: 
7 > 0.5. 

These inconsistencies can be resolved by looking at the 
scaling of the maximum vertical velocity in the simula- 
tions, where it is found that Re mnr ~ Ra o b for all three 
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simulated cases. Furthermore, this maximum is always 
found roughly midway up a sidewall in a strong persistent 
plume. Midway up a sidewall is also where the experi- 
mental measurements of velocity are taken. Therefore, 
this indicates that these experimental velocity measure- 
ments are not representative of the velocity as a whole, 
but only representative of the velocity at the walls that 
probably comes from the maximum possible vertical ve- 
locity within plumes. The reason AR = 4, Pr = 7 gives 
7 > 0.5 can be understood by noting that even at its 
highest Ra of 10 7 , the Re for this case is not turbulent 
and that in the published visualizations this case 
is dominated by laminar plumes. Therefore, the veloc- 
ity in individual plumes seems to obey ^ ma x = 0.5, but 
when the flow is turbulent the average exponent is closer 
to 7 = 0.46. This could represent an average between 
Imax = 0.5 and the exponent found for the Reynolds 
number dependence of the large-scale circulation, where 
7 = 0.43 fjj. 



Exp on 


Expt. 


AR = l,Pr = 


0.7 4,0.3 


4,7 


/3t 


0.309 


0.27±0.02 


0.31±0.02 


0.30±0.02 


Tmax 


0.50 


0.50±0.01 


0.50±0.01 


0.51±0.01 


s c 


0.145 


0.17±0.02 


0.14±0.02 


0.16±0.03 


7 


0.454 


0.46±0.02 


0.46±0.01 


0.52±0.03 


Re(Ra = 


10 7 ) 


434 


3000 


14 



TABLE I. Experimental [2] and numerical exponents. As- 
pect ratio and Prandtl number for the three numerical cases 
is given. The experimental value for velocity scaling is given 
as 7max- How the experimental value for 7 is gotten is ex- 
plained in the text. Errors for the simulations are based upon 
scatter of exponents taken between Ra. Re = udjv is taken 
at Ra = 10 7 for all three cases. 



It has been found that the dissipation at the wall ap- 
proaches a multiple of the mean dissipation across the 
box, contrary to the assumptions of the effect of shear 
upon the thermal boundary layer. Therefore, it does not 
seem that the dynamics in the boundary layer can be 
governed by shears, whether they be laminar or turbu- 
lent. Instead, it is suggested that plumes might domi- 
nate the dynamics and perhaps a new model of turbulent 
convection should be constructed based upon plume dy- 
namics. The basis for this suggestion is how the velocity 
scales. Evidence is presented that the experimentally ob- 
served scaling of the Reynolds number as Re ~ Ra 1 , with 
7 = 0.5 can be reproduced by the simulations only if this 
is taken to be the maximum vertical velocity in persis- 
tent plumes along the sidewalls. So one should replace 7 
by 7max in the experiments. Based upon the simulations, 
it is proposed that the simulated exponent based upon 
the average kinetic energy is 7 sa 0.46, which is shown 
to be consistent with the experimental measurements of 
the scaling of the temperature fluctuations and Nu. This 
suggests that it would be useful to have several determi- 
nations of the Reynolds numbers from a single experi- 
ment. That is, Reynolds numbers based upon velocity 
measurements taken near a sidewall, in the interior, and 
from the large-scale circulation. 
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